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Abstract: More than 25 years ago, Mueller Navelet jets were proposed as a decisive test 
of BFKL dynamics at hadron colliders. We here study this process at NLL BFKL accu- 
racy, taking into account NLL corrections to the Green's function and to the jet vertices. 
We present detailed predictions for various observables that can be measured at LHC in 
ongoing experiments like ATLAS or CMS at ^/s = 7 TeV: the cross-section, the azimuthal 
correlations and the angular distribution of these jets. For this purpose, we apply realistic 
kinematical cuts and binning, and study the dependence of our results with respect to 
several parameters. We then compare our results with those that can be obtained in a 
fixed order NLO treatment, and propose specific observables which could actually be used 
as a probe of BFKL dynamics. 
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1. Introduction 

One of the important longstanding theoretical questions raised by QCD is its beliaviour in 
the perturbative Regge hmit s S> —t. Based on theoretical grounds, one should identify 
and test suitable observables in order to test these peculiar dynamics. 

First, one should select processes in which the presence of a hard scale justifies the use of 
perturbative QCD. At high energy, QCD is a massless theory with vector bosons, and it has 
two kinds of infrared (IR) divergences, namely the soft and the collinear divergences. For 
sufficiently inclusive quantities, both kinds of divergences cancel. Still, they are responsible 
for large logarithms, which may compensate the smallness of the strong coupling. At 
leading order, the soft singularities manifest themselves as powers of as In s/\t\, resummed 
by the leading logarithmic (LL) BFKL Fomeron [1-4]. The collinear singularities are 
responsible for large logarithms of ratios of the transverse scales, which are resummed at 
leading logarithmic order (LLQ) by the DGLAP equation [5-8]. 

The Regge limit is expected to be governed by the soft perturbative dynamics of QCD, 
which we want to reveal, and not by its collinear dynamics. The key point is thus to select 
processes in which the hard collinear scales are of similar magnitude, in such a way that 
the difference between a fixed order calculation and a collinear resummed result should be 
tiny, while a BFKL type of resummation should modify the predictions dramatically. 

During the last 25 years, there have been many attempts to see manifestations of 
BFKL resummation effects. In inclusive DIS at HERA [9-11] or in total 7*7* cross-section 
at e+e^ colliders [12-18], the hard scale is the 7* virtuality. Exclusive processes have also 
been proposed and studied, either for heavy meson production {J/^, T), the hard scale 
being provided by the meson mass [19-22], or meson electroproduction at large t [21,23], 
for which HERA data seems to favour a BFKL picture [24, 25] . At future high energy and 
high luminosity colliders like ILC, processes like P P could be a realistic exclusive 

test of the hard Pomeron [26-31], with the planned detectors designed to cover the very 
forward region. 

Jets have been proposed as a powerful tool in order to study BFKL dynamics, like 
diffractive high energy double jet production [32-34] as well as central jets [35,36] in hadron- 
hadron collisions. In this paper, we focus on Mueller-Navelet jets [37]. This test of BFKL 
is based on the measure of two jets at large pr (hard scale), such that s S> pf- 3> Aq^jj, 
separated by a large rapidity Y , including possible activity between the two observed jets, 
as illustrated in figure 1. The idea is to consider two jets of similar pT in order to minimize 
the effect of collinear resummation. From a lowest order treatment it is clear that these 
two jets should be back-to-back, in the very forward and very backward regions. 

On the other hand, the expectation is that the large value of y = ln(s/p|,) should 
examplify the effect of BFKL dynamics, due to possible emission of gluons between them 
(thus the Pomeron contributes there at t = at the level of the cross-section), leading 
to enhanced terms which sum up as ^(a^y)" (LL), as^ia-gY)'^ (NLL) [38,39], etc., 
leading to a power-like rise for the cross-section. 
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Figure 1: kinematics 



Besides the cross-section a more exclusive observable within this process drew the 
attention, namely the azimuthal correlation between these jets [40,41]. The signal of a 
BFKL dynamics is a decorrelation of relative azimuthal angle between emitted jets when 
increasing Y . Indeed, while a fixed order calculation implies that the two jets would be 
emitted back-to-back, the fact that more and more (untagged) gluons can be emitted 
between them when increasing their relative rapidity should lead to a decorrelation of 
this relative azimuthal angle. Studies were made at LL [40-42], which overestimates this 
decorrelation by far. A better agreement with the data [43] could be obtained in the LL 
scenario using an event generator which takes into account in an exact way the energy- 
momentum conservation, which is a subleading effect in a pure BFKL approach [44]. On 
the other hand, the (kinematically) modified LL BFKL approach [45], again based on LL 
jet vertices, could also provide some better agreement with the data. 

At the same time, an exact fixed NLO (af) Monte Carlo calculation using the program 
JETRAD [46] lead to a too low estimate of the decorrelation, while the Monte Carlo NLO 
program HERWIG [47] was in perfect agreement with the data. It should be noted that 
this last treatment includes some Sudakov resummation effects, which might be important. 
The inclusion of such effects within a BFKL approach in an open problem which might be 
of interest for phenomenology. We leave this issue for further studies. 

Starting from first principles from the point of view of Regge and Quasi-multi-Regge 
kinematics, NLL [48, 49] and collinear resummed NLL [50] studies (with LL jet vertices) 
have been performed, improving the situation with respect to pure LL BFKL, but still 
leading to a much stronger decorrelation than the one seen by the data. 

In a previous work, we showed, based on a full NLL analysis [51], that contrarily to 
the expectation, the NLL corrections to the Green's function and to the jet vertices [52,53] 
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are of similar magnitude, based on a Mathematica code. We focused there on a center- 
of-mass energy ^/s =14 TeV, and considered jets with fixed values of transverse momenta 



In the present paper, we pursue this study and make detailed predictions for obser- 
vables to be extracted in the ongoing experiments ATLAS and CMS. Since experimental 
data are given in bins, it thus requires that we integrate kj j over a finite range. We also 
fix the ccnter-of-mass energy to be y/s =7 TcV. In the study of the azimuthal correlations, 
we extend the use of the collinear resummation method to non-zero conformal spins. We 
then study in detail the azimuthal distribution of jets, which is directly experimentally 
accessible. Finally, we make a detailed comparison (with the same set of parameters) of 
our predictions with the NLO fixed order results based on the code used in ref . [54] . 

Our numerical predictions are based on a new Fortran code which allows us to 
perform more detailed studies of the dependency on various parameters (PDFs, renormal- 
ization/factorization scale, choice of sq scale). To check the consistency of our results, a 
detailed comparison has been made, in the mixed case of NLL BFKL Green's function and 
LL jet vertices V^^\ with the previous studies of ref. [49,50]. 

2. Basic formulas for LL and NLL calculation 
2.1 Kinematics and general framework 

We consider two hadrons (in practice protons) which collide at a center-of-mass energy ^/s 
producing two very forward jets, whose transverse momenta are labeled by Euclidean two 
dimensional vectors kj,i and kj,2, and by their azimuthal angles and (f)j^2- The jet 
rapidities yj^i and yj^2 are related to the longitudinal momentum fractions of the jets via 

xj = -i^e^-^. The two partons produced by each of these two hadrons, which initiate the 
hard process, are treated in a collinear way. For large xj^i and xj^2i collinear factorization 
leads to a differential cross-section which reads 



proton, characterized by their longitudinal momentum fraction Xj. The hard process is 
then described using fcr-factorization. The logarithmically enhanced contributions are 
taken care of by convoluting, in transverse momentum space, the BFKL Green's function 
G with the two jet vertices, according to 



where the Mandelstam variable s = X1X2S refers to the hard subprocess. The jet vertices 
Va^b were calculated at NLL order in ref. [52,53]. Combining the PDFs with the jet vertices. 





d</'j,id0j,2 / d2kid2k2l4(-ki,a;i)G(ki,k2,s 



)Vb{k2,X2) , 

(2.2) 
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we can thus write 



where 



= J d(/>j,id(/.j,2y^d2kid2k2$(kj,i,Xj,i,-ki)G(ki,k2,l)$(kj,2,2;j,2,k2), (2.3) 

^(kj,i,xj,i,ki) = J dxif{xi)V{lLi,Xi). (2.4) 

In order to deal both with the cross-section and with the azimuthal decorrelation, it is 
convenient to define the coefficients 

Cm(|kj,i|, \kj^2\,Y) = J dyi dy2 5{yi +y2-Y) J d^j,i d^j,2 cos {micpj^i - (pj^2 - vr)) 

X y"d2kid2k2$(kj,i,xj,i,-ki)G(ki,k2,s)$(kj,2,xj,2,k2). (2.5) 

The differential cross-section then corresponds to Cq which reads 

[ dyidy2 5{yi + y2-Y)— = — — = Cq , (2.6) 

J d|kj,i|d|kj,2|dyj,idyj,2 d|kj,i| d|kj,2| cZF 

while the azimuthal decorrelation for fixed (|kj^i|, |kj^2|,^) is given by 

C 

{cos{m(p)) = (cos (m(0j,i - (pj,2 - ^r))) = 7^ • (2.7) 

2.2 LL order 

The jet vertex V at lowest order just implements the fact that the jet is made of a single 
parton, of the same nature as the collinear parton initiating the hard process. It reads [52, 
53]: 

Fi°)(k,x) =/.i°)(k)5f (k;x), /^i°)(k) =^%^, (2.8) 



Sf\k;x) =s(l-^) |kj|<5(2)(k - kj). (2.9) 



X 

In the definition of hi^\ Ca = Nc = 3 is to be used for initial gluon and Cf = {N^ — 
l)/{2Nc) = 4/3 for initial quark. 

In the LL approximation, the BFKL kernel, because of its conformal invariance, is 
diagonalized by the eigenfunctions 



E„,,(k,) = ^(kf)'^-^e-^^ (2.10) 



with an eigenvalue given by 



uJLLin, v) = asXo + , (2.11) 
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with as = Ncas/n and 

Xo(n, 7) = 2*(1) - * (7 + I) - * (l - 7 + I) , (2.12) 

where ^{x) = r'(a:)/r(a:). Using this basis for both the Green's function and the jet 
vertices, one thus obtains, introducing the arbitrary (at LL) scale sq, 



Cm = (4-35^,0) J dyi dy2d{yi+y2-Y)J di^ Cm,ui\^j,i\, xj^i)C^^^,{\kj^2\, xj. 
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, u!(m,iy) 



(2.13) 
where 

Cm,.(|kj|,xj) = y d(^jd2kdx/(x)F(k,x)£;m,.(k)cos(m<^j). (2.14) 
2.3 NLL order 

At NLL, the jet can be made of either a single or two partons. The coUinear singularities 
can be absorbed consistently in the renormalizcd PDFs, as was shown in refs. [52,53], for 
a given infrared-safe jet algorithm. These jet vertices read symbolically 

K(k,x) = yi°)(k,x) + a,FW(k,x). (2.15) 

The explicit form for the NLL vj^^^ are rather lengthy and are intimately dependent on 
the jet algorithm. They will not be reproduced here^. In our study we will use the cone 
algorithm with a size of -Rcone = 0.5. Note that other jet algorithms can be used, which do 
not affect significantly our obtained results. 

The main issue when dealing with NLL corrections is to treat the NLL BFKL kernel. 
The point here is to avoid dealing explicitly with two convolutions in transverse momentum 
space, between the jet vertices and the Green's function. In general this kind of convolution 
is very difficult to handle with for numerical evaluations. Instead we prefer to mimic the 
treatment used for LL studies and work in the (n, u) space. One is thus looking for a 
convenient basis in order to deal with the NLL BFKL kernel. The functions (2.10) cannot 
be used in principle, since conformal invariance is now broken. Anyway, the action of 
the NLL BFKL kernel on these LL eigenfunctions has been calculated in ref. [56], and it 
turns out that En^i, are still eigenfunctions in an extended sense, if one now promotes the 
eigenvalue to become an operator containing a derivative with respect to u [29,48,49] . When 
convoluting with jet vertices, this derivative acts on them, thus leading to a contribution 
to the eigenvalue which now depends on the jet vertices [29, 48, 49, 57] 



i^NLL{n, v) = asXQ ( I'^l, ^ + ) + «^ 



Xi ( \n\,^ + iv 



-Xo |?^|, o i -21n/x^-?— In 



2A^e"" V 2 ' ; i ^« du C„,,(|kj,2|,Xj,; 



(2.16) 



^They can be found in ref. [51], as extracted from refs. [52,53] after correcting a few misprints of ref. [52]. 
They have been recently reobtained in ref. [55]. 
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where 



3 Pq 



+ 



r (7 + ^) + < (1 - 7 + f ) - 2</>(n, 7) - 2</>(n, 1 - 7) 



TT^ cos(7r7) 



4sin^(7r7)(l - 27) 



3+1 + 



- 1 



Nf\ 2 + 37(1-7) 
N^J (3-27)(l + 27) 

7(1-7) 



iV3y 2(3-27)(l + 27) 



(2.17) 



with the constant 5 = (4 - tt^ + 5Po/Nc)/12. ({n) = Efeli is the Riemann zeta 
function while the function (f) reads 



(fc + n + l)-V'(fc + l) 



+ (-1)^+1 [;3'(A; + n + 1) + ;5'(A; + 1)] + 



V'(fe + i) -V(fe + " + i) 



A: + 7+f 



with 



V'' 



1 + 7 



(2.18) 



(2.19) 



At NLL accuracy, only the leading order vertex coefficients enter in the derivative term 
of (2.16), so that 



2 Q 

—2 In Ud — i— In 
av 



d^°)(|kj,i|,xj,i) 



= 21n 



(C£?)(|kj,2|,xj,2))* 

2.4 Strong coupling, renormalization scheme and PDFs at NLL 



(2.20) 



In this paper we will mainly use the MSTW 2008 PDFs [58]. We will make comparisons 
with several other sets of PDFs, as provided by the Les Houches Accord PDF Interface 
(LHAPDF) [59]. 

We use the two-loop strong coupling constant in the form 



bi InL 



with L = In/Ltlj/AqQj-), and 



60 = 



33 - 2Nf 



127r 



61 = 



153 - 19iV/ 



247r2 



(2.21) 



(2.22) 



In the following, or without argument is to be understood as a5.(/i|j) or as{fi^) 
respectively. The MSTW 2008 PDFs assume ^u/j and ^u^? to be equal. Therefore, we make 
the same identification everywhere in our analysis. The renormalization scale /ir is chosen 
to be i^R = A/|kj,i| • |kj,2|. 
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2.5 Choice of scale sq 



At NLL, one should also pay attention to the choice of scale sq. The choice sq = ^so,i *o,2 
with so,i = -^kj • which we adopt is natural, since it does not depend on the momenta ki^2 
to be integrated out. Besides, the dependence with respect to sq of the whole amplitude can 



be studied, when taking into account the fact that both the NLL BFKL Green's function 
and the vertex functions are so dependent. We refer to Sec. 3.2.2 of ref. [51] for a detailed 
discussion. 

2.6 Collinear improvement 

Several methods have been developed to improve the NLL BFKL Green's function for 
n = 0, by imposing compatibility with the DGLAP equation [5-8] in the collinear limit 
[60-63]. This is only required by the Green's function. Indeed, the collinear improvement 
deals with poles in the 7 plane (7 being a variable conjugated to transverse momentum 
in the Mellin transform). We have checked, based on a numerical study [51], that the 
jet vertices are free of 7 poles and thus do not call for any collinear improvement. In 
order to study the effect of such possible collinear improvement [60-63], a first attempt 
was performed in ref. [51], for n = 0, using the scheme 3 of ref. [60]. Focusing on n = is 
enough for the study of the cross-section. 

In view of the study of azimuthal correlation, a consistent treatment requires to take 
into account these collinear improvements also for n 7^ . This has been investigated for 
the NLL BFKL Green's function in refs. [49,50,57]. We take into account these effects in 
this paper, thus improving the study of our previous work [51]. 

3. Binning in |kj| 

3.1 Integration over |kj| 

The experimental binning imposes that the values of |kj^j| should be integrated in a given 
range. Each rapidity yi varies in the range ymin ^ Ui ^ Umax- In practice, we take ?/min = 
and ymax = 4.7 . The total relative rapidity Y = yi + y2 which is experimentally accessible 
varies between and 9.4. We will restrict ourselves to the region Y > A (see discussion at 
the end of this section). 

The phase space, at fixed |kj^i|, |kj^2|, Y, is defined as 



d(P. S) = d|kj,i| d|kj,2| dyi dy2 S{yi +y2-Y). 



(3.1) 



Correspondingly, the integration over the bin phase-space is defined as 




(3.2) 



For a given observable O , we thus define 




(3.3) 
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which equivalently can be written as 



''J max, 2 



Obi, 



J,i\ 



J d\kj^2\J dyi J <5(yi+y2-5^)0(|kj,i|,|kj,2Ui,y2) (3.4) 



^Jmin,l ^.7min,2 
^.7max,l ^Jmax,2 1 



2/min Vmin 



j (i|kj,i| j d\kj^2\ J dyi 9(ymin < - yi < ymax) Odkj,!], |kj,2|,yi,y - yi) . (3.5) 



The resulting cross-section 



da 
dY 



bin 



bin d|kj,i| d|kj,2| dyj,i dyj,2 



(3.6) 



is in practice numerically evaluated by sampling each i/i with a (ymax — 2/min)/10 = 0.47 
binning. It thus means that in eq. (3.5), the yi integration is replaced by a discrete sum, 
which is then multiplied by a 0.47 width. This cross-section is shown in figure 2, in the 
pure LL approximation as well as in the full NLL treatment. This figure shows that a very 
significant fraction (~ 80%) of the cross-section is obtained for A;jmax ~ 60 GeV. We will 
further discuss this in the next subsection in relation with energy-momentum conservation 
issues. 




160 



^Jmax [GeV] 



Figure 2: Growth of the cross-section with fcj, 

fcjmin,! = fc.7min,2 = 35 GeV. 



= fcjmaxa = fc,/max,2 , for F = 6.6 and fcjmin = 



As long as the jet vertices are treated in the LL approximation, the integration with 
respect to |kj^j| can be performed analytically. As a consistency check, we compare, in the 
Tevatron kinematics used in ref. [49], the integration with respect to [kj^jj with boundaries 
fejmin,! = 20 GeV, /cjmin,2 = 50 GeV and /cjmax.i = ^Jmax,2 = oo . Numerically, due to 
numerical instabilities when evaluating the Green's function for fixed values of [kj^jj at low 
Y, the comparison with data is expected to be rather poor for Y < 2asNc ' ^yP^^^^^^y 
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for y < 4 . A detailed study of this issue will be made elsewhere [64] . In the rest of this 
paper, we will restrict ourselves to the region Y > 4 . 

To perform this comparison, we use the observable (cos , which is the average of 
(cos if) on the experimental bin, that is defined here as (see eq. 2.7)) 



(cos(mv9))l,i 



"'Jmax,! J max, 2 

/ ,/ 

^Jmin.l ^' ./mill, 2 



d|kj,i| d|kj,2|C,i 



bin 



^Jmax,l 'i'Jmax,2 



I I d|kj,i|d|kj,2|Co 



^ Jmin.l ^ Jmin,2 



da 



d(P.S)(cos(m(/j)) — 

bin d|kj,i|d|kj^2|d7/j,id7/j,2 

^^^d{P. S) ^jj^^^^i d|kj^2| dyj^i dyj,2 



(3.7) 



This comparison is shown in figure 3, showing the consistency of our numerical results^. 



(cos if) 



bin 




pure LL analytical 

+ pure LL numerical 

NLL (Green's function) + LL (jet) analytical 

X NLL (Green's function) + LL (jet) numerical 



Y 



Figure 3: Comparison of (cos(p)^^^, either using the numerical integration over jkj^l, |kx2|, based 
on our Fortran code (crosses), or the analytical integration (solid) as shown in ref. [49], for the 
Tevatron kinematics, for both pure LL BFKL and mixed NLL BFKL treatments. 



3.2 Energy-momentum conservation issues 

It is well known that the BFKL equation does not preserve energy-momentum conservation. 
However, this violation is expected to be smaller at higher order in perturbation theory, 
i.e. when comparing NLL BFKL versus LL BFKL. In practice, one should thus avoid to 
use all the available collider energy. This means that one should satisfy the constraint 

yj,, < cosh-i ^ . (3.8) 

^This comparison focuses on the Green's function. It thus assumes that the PDFs are equal to 1, 
with xj^i = 1. The scales fj,F and y'so are taken to be y^kjmin, i kjmin, 2 , with fcjmin, i ~ 20 GeV and 
fcjn.in,2 = 50 GeV. 
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This implies that taking a lower kj provides a larger validity domain. This justifies in our 
opinion a strong experimental effort to extract low kj data. 

In practice, with only a lower cut on kj, one has to integrate over regions where the 
BFKL approach may not be valid anymore. For example, kj = 60 GeV leads to a constraint 
yj,i ^ 7.3 . For this reason it would be nice to have a measurement with bins including 
an upper cut on transverse momentum, fcjmin < fcj < ^Jmax- Since the cross-sections 
are expected to be large, we believe that the statistics should be large enough to allow 
for a narrow binning in kj, which should thus be mainly a detector issue"^. To conclude 
this section, we note that the kj integration reduces the Y domain between jets. This Y 
domain is also reduced by the Xj integration, which is weighted by PDFs, strongly peaked 
at small Xi . 



4. Results: symmetric configuration 

In this section, we consider a symmetric configuration as planned to be studied by the 
CMS collaboration. We thus consider bins with cuts 

35 GeV < |kj,i|,|kj,2| < 60 GeV, 

0< 2/1,2/2 <4.7. (4.1) 

We consider several kinds of scenarios, starting from a pure LL approximation up to full 
NLL and collinear improved NLL approximations. The convention for colors is the same 
in the whole paper: 

blue: pure LL result 

magenta: combination of LL vertices with pure NLL Green's function 

green: combination of LL vertices with collinear improved NLL Green's function 

brown: pure NLL result 

red: full NLL vertices with collinear improved NLL Green's function. 

(4.2) 

4.1 Cross-section 

We first consider the cross-section. The obtained results are displayed in figure 4. Note 
that the Monte Carlo integration leads to a precision of the order of 2% to 5%, which is 
too small to be seen in this figure. 

This result confirms the fact that NLL corrections to the jet vertices are huge, of the 
same order of magnitude as the NLL corrections to the Green's function. The full NLL 
result leads to a cross-section which is significantly smaller than the one based on LL 
vertices combined with the pure NLL Green's function. 

The curves obtained when combining the LL vertices with the pure NLL Green's 
function and when combining the LL vertices with the collinear improved NLL Green's 



^At CMS, a measurement with fcjmin = 35 GeV seems to be possible, while going down to 20 GeV 
would probably require a dedicated trigger. 
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function are almost indistinguishable. Similarly, the curves obtained when combining the 
NLL vertices with the pure NLL Green's function and when combining the NLL vertices 
with the collinear improved NLL Green's function are very close. 

In figure 5, in the pure NLL case, we display the uncertainties due to the changes of 
the various involved parameters. The first effect which we study is the variation of the 
scales So and /^i?, as shown in figure 5 (L). The large uncertainty at very low Y is related 
to the specific instabilities of NLL Green's function mentioned at the end of section 3.1, 
while the large Y uncertainty is related to kinematical boundary effect (the cross-section 
almost vanishes). 

The second effect, due to the dependency on the set of PDFs, is shown in figure 5 (R). 
Note that we only display the pure NLL case, although the trend is similar for other 
scenarios. 

Both of these effects are much smaller than the changes due to the NLL corrections to 
the jet vertices. 




Y 



4 5 6 7 8 9 



Figure 4: Differential cross-section as a function of the jet rapidity separation Y, integrated over 
bins 35 GeV < |kj_i|, |kj^2| < 60 GeV and < yi, y2 < 4.7, for the 5 scenarios described in the text, 
see (4.2). 

4.2 Azimuthal correlations 

We now consider the azimuthal correlations. The obtained results are displayed in figure 6, 
again using the color conventions (4.2). Note that the Monte Carlo integration leads to a 
precision of the order of a few % when using the NLL vertices, the numerical uncertainty 
being negligible in the case of LL jet vertices. We do not show it on this figure. 

Comparing on one hand the pure LL scenario with the mixed LL vertex combined 
with the NLL Green's function, and the mixed LL vertex combined with the NLL Green's 
function with the full NLL treatment on the other hand, we see that the correction due to 
the jet vertex produces the largest correction. On the same plot, one can see the effect of 
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Figure 5: Left: Relative variation of the cross-section when varying ^/sq and with a factor 2. 
Right: Relative variation of cross-section with respect to MSTW PDFs due to the replacement by 
other sets of PDFs, as indicated. 




Figure 6: The bin averaged {cosip)^^^^ as a function of the jet rapidity separation Y, integrated 
over bins 35 GeV < |kj^i|, |kj_2| < 60 GeV and < yi, 2/2 < 4.7, for the 5 scenarios described in the 
text, see (4.2). 



collinear improvement. When including this effect for the whole set of the conformal spins 
n, we obtain very close results for the pure NLL and the collinear improved NLL approaches. 
This can be compared with the mixed LL jet combined with either a pure NLL Green's 
function or collinear improved NLL Green's function scenarios: the resulting modification 
is much smaller at NLL. This is in our opinion a sign of convergence of the perturbative 
series. One can see on figure 7 this effect already at the level of fixed |kj^i| = |kj.2| = 35 
GeV, when passing from n = collinear improved NLL Green's function (left) to all n 
collinear improved NLL Green's function (right)^. Note that this resummation affects only 

*We plot these curves for ^/s — 7 TeV, with the same choice of parameters for yt as in ref. [51], i.e. 
3 < 2/i < 5 and thus 6 < F < 10. 
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the green and the red curves in figure 7, which are thus the only modified ones when 
passing from figure 7 (L) to figure 7 (R). Furthermore, we see that including the collinear 
resummation for all n does not lead anymore to (cos (p) potentially above 1 (except in the 
very large Y domain, due to the kinematical boundary effect). 

We have made a similar check, in the three scenarios with LL jet vertices, between 
the results obtained with our code and the one used in ref. [50]. These results are in very 
good agreement, taking into account the slightly different numerical treatments and the 
fact that we use scheme 3 of ref. [60] while the code of ref. [50] is based on scheme 4. 



(cos ip) 



(cos if) 




Y 




Figure 7: Left: collinear resummation for n — Q\ right: collinear resummation for all n. 



The predictions of figure 6 show that, contrarily to the natural expectation, the inclu- 
sion of mini-jets between the two tagged jets, when performed at full NLL, does not break 
the very high correlation between these two jets. Moreover, the obtained decorrelation 
effect is very flat with respect to Y . This is an effect which takes origin from the NLL 
corrections to vertices. 

We now study the stability of this result with respect to changes of parameters, within 
the pure NLL approximation. The variation due to change of the scales sq and /x^r is 
shown in figure 8 (L). The effect is sizeable, but does not change the conclusion that the 
decorrelation remains much smaller than in the pure LL or mixed LL-l-NLL approaches. 
The second effect, due to the dependency on the set of PDFs, is shown in figure 8 (R). 
This dependency is very weak, much smaller than for the cross-section, see figure 5 (R). 
Note that this PDF dependency does not exist when using LL jet vertices. 

Let us now consider the observable (cos(2(^))jjjj^. The results based on the 5 approaches 
(4.2) are displayed in figure 9 (L). Similar conclusions as for (cos ^)]^^^ can be drawn. Indeed, 
an even more dramatic effect due to the NLL corrections to the jet vertices is observed. 
On the other hand, the difference between the pure NLL and the collinear improved NLL 
treatments is very small, this time of the same order of magnitude as the one observed 
between a mixed LL jet with pure NLL Green's function and the mixed LL jet with 
collinear improved NLL Green's function approaches. 

Again, when including full NLL corrections, the decorrelation effect is rather small, 
and much smaller than the one obtained in non full NLL treatments, and the dependency 
with respect to Y becomes much more fiattish. 
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Figure 8: Left: Variation of {cos ip) t^-^^ when varying and witlr a factor 2. Riglrt: Relative 
variation of {cofnp)-^^^^ with respect to MSTW PDFs due to the replacement by other sets of PDFs, 
as indicated. 



In figure 9 (R), we show the dependency of our full NLL prediction with respect to 
changes of the scales sq and ^p- The variation due to sq changes is very small, almost 
negligible in comparison with the same dependency for (cos '/')bin (^^^ figure 8 (L)). Besides, 
the dependency with respect to fip remains sizeable, although a bit smaller in absolute 
magnitude than in {cos if) ^^^^ (but comparable in relative magnitude). 

Again, these dependencies do not change the conclusion that the decorrelation remains 
much smaller than in the pure LL or mixed LL+NLL approaches. The second effect, due to 
the dependency on the set of PDFs, is very weak, similar to the one shown in figure 8 (R), 
and will not be displayed here. 




Figure 9: Left: The bin averaged {cos{2(p)) ^-^^^ as a function of the jet rapidity separation Y, 
integrated over bins 35 GeV < |kj_i|, |kj_2| < 60GeV and < yi, y2 < 4.7, for the 5 scenarios 
described in the text, see (4.2). Right: Variation of {cos{2(p)) ^^^^ when varying y/so and with a 
factor 2. 

It turns out that this remaining dependency with respect to ^i;' is much reduced when 
considering the observable {cos{2ip))^-^^/ {cos if)^^^ . This observable is shown in figure 10. 
In figure 10 (L) we display our prediction based on the 5 approaches (4.2), while in fig- 
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ure 10 (R) we show the and dependency. The difference between the full NLL 
prediction (either collinearly improved or not) and the non-full NLL ones is sizeable for 
Y > 6 , and the figure 10 (R) explicitly shows that this remains valid when taking into 
account ^/sq and f^p dependencies. 
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Figure 10: Left: The bin averaged {cos{2ip))y^^^/ {cosip)^-^^ as a function of the jet rapidity sep- 
aration Y, integrated over bins 35 GeV < |k,7^i|, |kj2| < 60GeV and < yi, y2 < 4.7, for the 5 
scenarios described in the text, see (4.2). Right: Variation of {cos{2ip))^^^/ {cos ip)^^^ when varying 
^/sq and (ip with a factor 2. 



The extraction of higher harmonics can be as well experimentally performed. We 
show in figure 11 (L) our predictions for (cos(3v9))jjjj^ based on the 5 different treatments 
(4.2), and the corresponding sensitivity with respect to and fip in figure 11 (R). In 
comparison with {cos{2(p)) , again the effect of NLL corrections in jet vertices is very 
important, leading to a much smaller decor relation. The dependency with respect to sq is 
similarly small, while the fip dependency is still sizeable. It is smaller in absolute magnitude 
than the one for {cos{2ip)) , although comparable in relative magnitude. 
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Figure 11: Left: The bin averaged (cos(3(^))jjjj^ as a function of the jet rapidity separation Y, 
integrated over bins 35 GeV < |k,/^i|, |kj,2| < 60GeV and < yi, y2 < 4.7, for the 5 scenarios 
described in the text, see (4.2). Right: Variation of (cos(3i^))^^jjj when varying ^/so and with a 
factor 2. 



15 



In figure 12, we show predictions for the observable {cos{3ip))^^^/ {cos{2if))^^^ , which 
is less sensitive to changes of factorization scale /ip- In figure 12 (L) we see that the 
differences between the approaches (4.2) are not sizeable. We see in figure 12 (R) that 
this remains true when taking into account the ^/sq and fip dependency. 
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Figure 12: Left: The bin averaged {cos{3ip))-^-^^/ {cos{2ip))^^^^ as a function of the jet rapidity 
separation Y, integrated over bins 35GeV < |kj^i|, |kj_2| < 60GeV and < yi, y2 < 4.7, for the 
5 scenarios described in the text, see (4.2). Right: Variation of {cos{3^p))^^i^^/ {cos{2ip))^.^^ when 
varying y/so and hf with a factor 2. 



4.3 Azimuthal distribution 

In practice, the observable which is most directly accessible in experiments is the azimuthal 
distribution of the two jets, defined as 

1 dcj 1 f °° 1 

^ — |l + 2^^ cos (nv9) (cos (nv9))| . (4.3) 

It is shown in figure 13 for LL, mixed jet LL with NLL Green's function, and mixed jet LL 
with collinear improved NLL Green's function, and in figure 14 for pure NLL and collinear 
improved NLL approaches. 

The figure 13 shows that the inclusion of NLL corrections to the Green's function 
leads to a smaller decorrelation compared to a pure LL treatment. Comparing figure 13 
with figure 14, we see that the NLL corrections to the jet vertices lead to an even larger 
correlation, at fixed Y. When increasing Y , we can see on these plots that the decorrelation 
effect is slower, as expected in BFKL picture. 

We now integrate also over Y in the range 6 < Y < 9.4 in addition to |k |kj^2|, yii 2/2- 
The resulting azimuthal distribution is shown in figure 15 for LL, mixed jet LL with NLL 
Green's function, and mixed jet LL with collinear improved NLL Green's function, and in 
figure 16 for the two full NLL scenarios of (4.2). In the same plots are shown the depen- 
dency with respect to sq and . We see that the pure LL approach is quite dependent 
on the scales and hf, whereas a mixed treatment using LL vertices with NLL Green's 
function shows a smaller dependency on np- The full NLL approaches are much more 
stable with respect to ^/sq, while still fip dependent. 



16 



Ida 

a dip 



1.4 

1.2 
1 

0.8 
0.6 
0.4 
0.2 




\da_ 

(7 dip 




~i r 

Y = 6.2 

Y = T.2 

Y = 8.2 




Figure 13: The azimutlial distribution for the 3 first scenarios of (4.2), for 3 values of Y. 
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Figure 14: The azimuthal distribution for the 2 full NLL scenarios of (4.2), for 3 values of 
Y. 



5. Results: asymmetric configuration 

An asymmetric configuration with very different kj^i^^ i and Arjmin. 2 allows us to compare 
our predictions with the ones obtained by fixed order NLO approaches, since it is known 
that symmetric configurations lead to unstable predictions in fixed order calculation [65,66]. 
Here we compare our predictions with the results obtained by the Duet code [54]. Below 
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Figure 15: The azimuthal distribution integrated over bins 35 GeV < |kj,i|, |k,7,2| < 60GeV, 
< yi, 2/2 < 4.7 and 6 < y < 9.4, for the 3 first scenarios of (4.2), including a variation of and 
with a factor 2 with respect to the central values. 
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Figure 16: The azimuthal distribution integrated over bins 35 GeV < |kj i|, |kj 2| < 60 GeV, 
< j/2 < 4.7 and 6 < F < 9.4, for the 2 full NLL scenarios of (4.2), including a variation of i^/sq 
and with a factor 2 with respect to the central values. 



we show the same observables as the ones which we considered in section 4 for the symmetric 
configuration, now supplemented by a comparison with the DiJET predictions, for which 
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we include a scale uncertainty on np of a factor 2, for every plot . We consider bins with 
cuts 

35GeV< |kj,i|,|kj,2| <60GeV, 
50GeV < Max(|kj,i|,|kj,2|), 

0< yi, 2/2 <4.7. (5.1) 

Since the cross-section is dominated by minimal allowed values of |kj^i|, |kj^2| ; such a choice 
of binning reduces the domain where |kj^i| and |kj^2| are very close to each other, for which 
unstable results at fixed order may be a source of worry. 

As the behaviour of the BFKL results is very similar to the one for a symmetric 
configuration, we will mainly focus in this section on the comparison with DiJET. 

5.1 Cross-section 

On the figure 17, we show the cross-section. This figure shows surprising results: the fixed 
order results are above BFKL predictions, contrarily to the expectation. This difference 




Figure 17: Differential cross-section as a function of the jet rapidity separation F, using cuts 
defined in (5.1). The different curves correspond to the 5 scenarios (4.2). The dots correspond to 
the predictions of the Duet code. 

remains valid after including the effect of possible variations of the parameters sq and // p , 
as can be seen in figure 18. 

5.2 Azimuthal correlations 

We now consider the azimuthal correlation, for which we again compare predictions based 
on 5 kinds of BFKL scenarios with a fixed order NLO prediction. The predictions for 

^Note that the results obtained with Duet use a scale fj,F = ^LLilililii^ ^ which is very close in the 
domain we consider to the value \/|kj.i| ■ |k,/,2| we use in our BFKL calculation 
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Figure 18: Differential cross-section as a function of the jet rapidity separation Y, using cuts 
defined in (5.1), in the full NLL approximation. The dots correspond to the predictions of the 
Duet code. We show the effect of a variation of ^/sq and ij,f with a factor 2 with respect to the 
central values. 



(cos(/3)j,j^ are displayed in figure 19 (L). The two full NLL BFKL predictions (pure and 
collinearly improved) are noticeably above the fixed order NLO prediction. However, the 
figure 19 (R) shows that the uncertainties with respect to sq and fip are quite significant, 
and do not allow to distinguish between the full NLL BFKL predictions and the fixed order 
NLO one. 
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Figure 19: Left: The bin averaged {cosip)-^-^^ as a function of the jet rapidity separation F, using 
cuts defined in (5.1), for the 5 scenarios described in the text, see (4.2). Right: Variation of 
(cos </?)i5jjj when varying and fxp with a factor 2. The dots correspond to the predictions of the 
Duet code. 



We then display the predictions for {cos{2(p))^-^j^ in figure 20 (L). The two full NLL 
BFKL predictions (pure and collinearly improved) are now a bit below the fixed order NLO 
prediction, and again, the figure 20 (R) shows that the uncertainties with respect to sq 
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and fip are quite significant, and do not allow to distinguish between the full NLL BFKL 
predictions and the fixed order NLO one with this observable. 




Y 



Figure 20: Left: The bin averaged {cos{2ip)) ^^^^ as a function of the jet rapidity separation Y, 
using cuts defined in (5.1), for the 5 scenarios described in the text, see (4.2). Right: Variation of 
(cos(2(/j))jjjj^ when varying ^/so and /ij? with a factor 2. The dots correspond to the predictions of 
the Duet code. 



Let us now study the observable 



(cos 'fi)hin 



The fact that in the full NLL BFKL 



predictions, {cos(p)^:^^ (resp. {cos{2(p)) y^^^) are above (resp. below) the fixed NLO order 
predictions now leads to a ratio which is very significantly, in the full NLL BFKL approxi- 
mation, under the NLO fixed order one, as can be seen from figure 21 (L). This difference 
is not washed out when including the uncertainties due to sq and variations, as shown 
in figure 21 (R). We want to emphasize the fact that this is valid in particular in the 
region y ~ 6, for which according to the discussion of section 3.2, the corrections due to 
energy- momentum conservation are not expected to be very significant. 
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Figure 21: Left: The bin averaged (cos(2(y9))|^jjy (cos 1^9)^^;^^ as a function of the jet rapidity sep- 
aration Y, using cuts defined in (5.1), for the 5 scenarios described in the text, see (4.2). Right: 
/(cos when varying 



Variation of (cos(2iy9))|^ 
to the predictions of the Duet code. 



So and fip with a factor 2. The dots correspond 
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We finally consider the observable (cos(3(/?))|-,;^ in figure 22 (L). The two full NLL 
BFKL predictions (pure and collinearly improved) are now significantly below the fixed 
order NLO prediction. The figure 22 (R) shows that the uncertainties with respect to sq 
and fip are quite significant, and can marginally alter the possibility of distinguishing the 
two types of scenarios, mainly due to the up uncertainty. Considering the ratio ^'^°'^(^'^))bin 

(cos (2>^)) bin 

we observe that this latter observable is much more favorable, as can be seen from figure 22. 
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Figure 22: Left: The bin averaged {cos{3(p)) ^^^^ as a function of the jet rapidity separation Y, 
using cuts defined in (5.1), for the 5 scenarios described in the text, see (4.2). Right: Variation of 
(cos(3(/3)){^jjj when varying and /li? with a factor 2. The dots correspond to the predictions of 
the Duet code. 



(cos(3y))t^i,^ 



1 1 


1 1 1 

I I I I ' ' : 


1 1 


1 1 1 



Y 



1 1 

^^-^ • • 


} 


} 




1 

I - 


pure NLL 










/^F A'f/2 










I-IF 2/iF 










~ y/SQ/2 










^/s(i -!■ 










• fixed order NLO 
1 1 








1 



Y 



Figure 23: Left: The bin averaged (cos(3i^))^jj^/(cos(2iy9))jjjj^ as a function of the jet rapidity sep- 
aration Y, using cuts defined in (5.1), for the 5 scenarios described in the text, see (4.2). Right: 
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We do not study here the azimuthal distribution as we have done in the symmetric case 
of section 4.3. Indeed it is not possible to confront our BFKL predictions with the fixed 
order NLO predictions of the Duet code, due to instabilities when evaluating the higher 
harmonics, which are necessary for a precise study of ip distribution, using the Duet code. 
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6. Limit of small-R cone 



A detailed study, based on the work of ref. [67] where the jet vertices were computed in an 
approximated small R treatment, shows that the difference between an exact treatment, 
as used in the present work, and that small R approximation is small. This is illustrated 
in figure 24 for R = 0.3 and in figure 25 for R = 0.5 , and shows explicitly the consistency 
of the two approaches. This small R limit has been used in ref. [68] for phenomenological 
studies^. In this paper, it is stated that sizeable differences are obtained when comparing 
with the results we got in ref. [51], with the same set of parameters. We believe that this is 
mainly due to the way NNLL corrections, which are beyond the precision of both studies, 
are treated. Indeed, when convoluting jet vertices with the Green's function, there is a 
freedom to neglect terms of magnitude a^Y . A close inspection on the way both papers 
deal with such contributions shows that in ref. [68], eq. (42), terms involving the product 
of NLL corrections in both vertices are explicitly neglected, while they are kept in ref. [51] 
and in the present study. 
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Figure 24: Relative difference between an exact treatment of the cone size and the small cone 
approximation for R — 0.3 and |kj_i| = |kj2| = 35 GeV , in the two full NLL BFKL scenarios. 
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Figure 25: Relative difference between an exact 
approximation for R = 0.5 and |kj.i| = |k,7,2| 35 
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treatment of the cone size and the small cone 
GeV , in the two full NLL BFKL scenarios. 



®We thanfc the authors of ref. [68] for pointing out that the vaiues for Ci and C2 given in ref. [51] should 
be multiplied by a factor of 2 to get the proper normalization. 
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7. Conclusions 



In this paper, we have continued our NLL analysis of Mueller-Navelet jets, at the par- 
tonic level, concentrating on the kincmatical conditions of ongoing experiments ATLAS 
and CMS at LHC. We have made a detailed study of the azimuthal distributions for the 
first time at full NLL BFKL accuracy. Finally, wc confronted our predictions with the 
fixed order predictions based on the Duet code. Our analysis does not take into account 
hadronization effects. We did not estimate the importance of potentially competing pro- 
duction mechanisms involving multiparton interactions. The evaluation of the importance 
of these effects is left for future analysis. 

The predictions of the present study confirm the main result of ref. [51] that the effect 
of NLL corrections to jets vertices is dramatically large, similar in magnitude to the one 
due to the NLL Green's function corrections. 

We have investigated the stability of our predictions with respect to changes of fac- 
torization scale fip, of scale sq and of sets of PDFs. For the cross-section, in comparison 
with scenarios using LL jet vertices, the predictions arc much more stable with respect to 
variation of and sq, and of similar small order of magnitude for PDFs variations. Our 
full NLL BFKL predictions are surprisingly sizeably below the fixed order NLO prediction. 

For the decorrelation effect, the full NLL BFKL predictions and fixed order NLO one 
are very close for (cos 99) and (cos2(/?). They are very fiat in rapidity Y, but still rather 
dependent on sq, and specially on , while weakly dependent on PDFs. 

We have taken into account the effect of collinearly improved NLL BFKL Green's 
function. Including these effects for non zero conformal spins n has an important impact 
on our predictions for azimuthal decorrelation. It leads to results which are very close to 
the pure NLL BFKL treatment. 

The angular (p distribution which we predict at full NLL BFKL is very strongly peaked 
around and does not evolve strongly with respect to Y . This prediction significantly differs 
from the ones based on LL jet vertices, and is stable when changing fip , sq and PDFs. 

Finally, we have shown that for the ratios (cos 299)/(cos and (cos 3(/?)/(cos 2(/?) the 
differences between NLL BFKL and fixed order NLO are sizeable, and stable with respect 
to scale choices. 

To conclude, our analysis suggests that the ratios of harmonics are most suitable 
observables to distinguish between full NLL BFKL predictions and fixed order NLO ones. 
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